Properties of dense partially random graphs 
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We study the properties of random graphs where for each vertex a neighbourhood has been previously 
defined. The probability of an edge joining two vertices depends on whether the vertices are 
neighbours or not, as happens in Small World Graphs (SWGs). But we consider the case where 
the average degree of each node is of order of the size of the graph (unlike SWGs, which are 
sparse) . This allows us to calculate the mean distance and clustering, that are qualitatively similar 
(although not in such a dramatic scale range) to the case of SWGs. We also obtain analytically the 
distribution of eigenvalues of the corresponding adjacency matrices. This distribution is discrete 
for large eigenvalues and continuous for small eigenvalues. The continuous part of the distribution 
follows a semicircle law, whose width is proportional to the "disorder" of the graph, whereas the 
discrete part is simply a rescaling of the spectrum of the substrate. We apply our results to the 
calculation of the mixing rate and the synchronizability threshold. 
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I. INTRODUCTION 

Many natural and artificial systems are composed of a 
large number of identical agents that interact. Examples 
of this kind of systems are numerous and widespread: 
collaborators in physics, neural networks, computer pro- 
grams (considered as a system of interacting subrou- 
tines), the World Wide Web, movie actors, etc. 

The interactions however need not be identical and can 
vary in strength and range. A simplified analysis of such 
a system can take into account only the pattern of inter- 
actions, abstracting everything else. What is left is usu- 
ally called a network and is mathematically represented 
by a graph (along this article both terms will be used 
interchangeably. ) 

A graph is composed of vertices (associated to the 
agents) connected by edges. Agents only interact with 
other agents if there is an edge joining the correspond- 
ing vertices, and the strength of the interaction is given 
by the weight of the edge. The interactions between 
graph theory and physical science, particularly physics, 
has been very fruitful Jj . 

In 1959 Erdos and Renyi started a whole new 
branch of Graph Theory by creating (and extensively 
studying) the concept of random graphs. These are 
graphs where each edge has a defined probability of being 
present, and this probability is independent of all other 
edges. They seem particularly well suited for the study 
of systems for which there is little information about the 
range of the interactions. In these cases it seems natu- 
ral to assign independent and equal probabilities to the 
different connections. 

The average distance of a graph is defined as the av- 
erage of the length (i. e. the number of vertices) of 
the smallest path joining two edges. The clustering co- 
efficient gives the average number of connections present 
between neighbors of a vertex, divided by the number 
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of possible connections within the neighborhood of the 
vertex. In their groundbreaking article of 1998, Watts 
and Strogatz Q (hereafter WS) showed that most real 
networks display a very short average distance combined 
with large values of the clustering coefficient. But ran- 
dom graphs do not fulfill these requirements because, 
even though the average distance is short, the clustering 
coefficient turns out to be rather small. Thus they pro- 
posed a new graph model, the small world graph (here- 
after SWG). To build a SWG one starts with a regu- 
lar and sparse substrate graph (with large clustering but 
also with large average distance), and then rewires some 
edges with a probability p, thus creating shortcuts. It 
was shown that a small number of shortcuts is enough to 
significantly lower the average distance while leaving the 
clustering coefficient almost unchanged. 

Since their proposal, the properties of SWGs have been 
intensively studied, with numerical as well as analyti- 
cal methods. Even though some numerical analysis have 
been performed Q , one of the properties that still resists 
analytical treatment is the spectrum of the adjacency 
matrix of SWGs. 

The constraint of sparsity is justified by the fact that 
many networks in nature display this characteristic. But 
one can also find systems where the pattern of connec- 
tions of every node spans a significant portion of the 
whole graph (i.e. the degree of the nodes is of the same 
order as the size of the network). Some examples of this 
are the network of train routes in India 's'l, the full re- 
action graphs of the metabolic network of E. Coli |^ Q 
and the network of the interacting units of the computer 
program Mozilla P| . 

In this article we study the properties of graphs for 
which the sparsity constraint has been dropped. This 
allows us to calculate analytically the average distance 
and the clustering coefficient, as well as the whole eigen- 
value distribution of the corresponding adjacency matri- 
ces. The spectra of these matrices can be considered as 
limiting cases of the ones of sparse SWNs. 

The values of the different properties calculated in this 
article are only averages over a certain family of graphs 
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(defined in the next section). Nevertheless, simulations 
support the idea that the properties of almost every 
graph of this family should tend, in probability, to the 
average values found. 

In Section we define the model and relate it to the 
Watts-Strogatz prescription. In section Hill we calculate 
the average distance and the clustering coefficient. In 
section IIVI we obtain the spectrum of the corresponding 
adjacency matrices and compare it to the spectrum of 
sparse graphs. A couple of applications of the results 
obtained are presented in Section In section IVTI some 
conclusions are drawn. 



II. THE MODEL 

A graph is a pair of sets {V,E), where V has N 
elements (or vertices) and E C V"^ has M elements 
(or edges). Two vertices vi and V2 are connected if 

{vi,v2) e E. 

WS Q proposed a graph model capable of interpolat- 
ing between order and randomness. The graphs in this 
model are built by taking a substrate, which is a graph 
displaying some regularity, and randomly rewiring some 
of its edges, by keeping fixed one end of some edges and 
redirecting the other end to a different vertex at random, 
but following some rule. Let us consider, for example, 
the unidimensional case, where each vertex in a ring is 
connected to fc/2 vertices to the left, and k/2 vertices to 
the right. The process begins by traveling clockwise on 
the ring, and for each vertex rewiring, with probability 
p, the connection that joins it with its first neighbor to 
the right. Once the circle is completed, a new round is 
made where now the connection rewired is the one to 
the second neighbor to the right. The process ends af- 
ter fc/2 rounds (because k/2 is the number of neighbors 
to the right). A similar process can be implemented for 
a higher dimensional substrate. The graph obtained is 
called a SWG. 

In our model of dense partially random graphs (here- 
after DPRGs), we only consider hypercubic lattices as 
substrates, with the hope that for other substrates things 
will not be very different, as is the case for sparse 
SWG ia]. In these lattices V C R''^ and each dimen- 
sion has a different connectivity parameter fc^, which 
means that each node is connected to a hypercube of 
k = Y[i=i fc' — 1 nodes. This defines a fc-neighborhood 
for each vertex. As opposed to the usual constraint of 
considering sparse networks, here we are concerned only 
with dense graphs, i. e. graphs where each vertex is con- 
nected to fc = 0{N) other vertices. Notice that, to have 
the same number of neighbors for all nodes, we are con- 
sidering periodical boundary conditions for the lattices. 
Thus, the one dimensional lattice is formed by points on 
a ring, the two dimensional lattice by points on a torus, 
etc. 

We randomize the graph in the following way: each 
edge in every fc-neighborhood is deleted with a probabil- 



ity 1 — pi and vertices that do not belong to the same 
fc-neighborhood are joined with probability p2. This is 
equivalent to saying that, on an empty graph, each ver- 
tex is joined by a short link with probability pi to every 
vertex in its fc-neighborhood and with probability p2 to 
those outside it, by a long link or shortcut. A graph 
generated with this prescription is called Gp-^p^ (7) , with 
7 = k/N. The family of all such graphs, for pi, p2 and 7 
fixed is called Gpip2{l)- With pi = 1 and P2 = one ob- 
tains the ordered substrate, whereas for pi = P2 = k/N 
one obtains a random graph 10]. 

Notice that each vertex has a different number of con- 
nections, whose average number is fcpi -|- {N — 1 — k)p2 ~ 
N {-fpi -I- (1 — 7)^2)- But in the limit treated in this arti- 
cle, of large values of N, the deviations from this average 
value are exponentially small. To interpolate with only 
one parameter between a fully ordered and a fully ran- 
dom graph, we fix this average value by requiring that 
7Pi -I- (1 — 7)p2 = 7 (along the paper, though, results 
will be displayed for general pi and p2, unless otherwise 
stated). Notice that in this case the average number 
of connections of each vertex is fc. Choosing pi as the 
preferred parameter, a graph satisfying the relationship 
shown above is called Gp^ (7) and the family of graphs 
with pi and 7 fixed is called Gpiij)- As mentioned, for 
Pi — J the family Qp-^ (|^3|) of completely random graphs 
is obtained. 

In our model the number of shortcuts is P2N, for large 
N. To be able to compare the graphs in Gpii'j) with 
SWGs, it is necessary to know the number of shortcuts 
of the graphs generated with the WS prescription. In the 
sparse case, it is known that this number is ~ pkN/2. 
This comes from the fact that the sparsity of the net- 
work ensures that when a link is selected, the probability 
that it will be rewired to a vertex inside the same fc- 
neighbourhood is ~ 1/N. But if the network is dense, 
this probability becomes nonvanishing. The following 
procedure provides a good approximation to the num- 
ber of shortcuts. 

Instead of disconnect and rewire the links sequentially 
in each round, let us assume that in each round pN short 
edges are deleted at once, and then pN random edges are 
added, which can be short or long. The process begins 
with only the S = kN/2 short edges present in the graph. 
We call St the number of short edges after round t and 
Lt the number of long edges. Notice that after deleting 
pN short edges, the probability that one of the random 
edges added is short is the quotient between the number 
of available short edges and the number of total available 
edges: {pN + S - St)/{L + S - St - U + pN) = (pN + 
S - St)/{L + pN), where L = N{N - l)/2 - S. We 
are using the fact that the number of edges is conserved 
{St + Lt — S). Using this, the evolution equations for 
large enough N can be obtained: 
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FIG. 1: Fraction of shortcuts as a function of the rewiring 
probabihty p. The symbols are averages taken on 30 matrices 
with A'^ = 6001. The lines are the theoretical predictions. 
Error bars are smaller than the symbols. 



The second summand in Eq. ^ corresponds to the edges 
deleted and the third to the fraction of edges added that 
are short. This formula can be iterated to obtain: 



S~L + L 



L 



L + pN 



t+i 



(3) 



After t = k/2 rounds, and in the hmit of large k and 
N such that 7 = k/N, we obtain: 



needed to reach any other vertex. This in turn imphes 
that the average distance in such a graph tends to 2 — p. 
For a general graph it is known that 2— p, with p = 
2M/N{N + 1) is a lower bound to the average distance. 
It is interesting to notice that this bound is achieved by 
some graphs (stars, for example). Random graphs, on 
the other hand, only achieve it in the limit of infinite 
N. In our model this implies that the average distance 
satisfies d>2 — pij — p2{l — 7). 

The graphs in Gpip2 (7) can be generated as the union of 
two random graphs. One of them is simply a completely 
random graph with edge probability p2 , called Gp^ ■ For 
the other, the probability that an edge is present is 
if it is a long edge and p — ^^Zp^ if it is a short edge. 
This union results in a graph equivalent to the one ob- 
tained by adding a number Na of edges to Gp^ , such that 
Na = (1 — P2)P1- Now every one of these additional 
edges has the effect of decreasing the average distance 
of the graph by at least M~^. Thus, averaging over all 
graphs in Qp-i^p2{l), the mean average distance must sat- 
isfy d < 2 - p2 -7P(1 -P2) = 2 -pi7 -p2(l - 7)- This, 
together with the upper bound, implies that, in mean, 
d = 2 — P17 — P2(l — 7) for large N . The fact that this 
value is a lower bound for general graphs ensures that, 
for large values of iV, almost every graph has an average 
distance equal to the mean. 

Notice that for graphs in Qp-^ (7), in the infinite TV limit, 
the average distance for pi < 1 is only a function of 7: 
d = 2 — 7. For pi — 1 and p2 = 0, the graph is circulant, 
and it is not difficult to see that the average length is 
(1 -I- 7)727 9]. In figure|21it can be seen that for finite 
but large values of N , the average distance falls very 
rapidly with the number of shortcuts. 



pr = ^./5^1-i^(l-exp(-^)) (4) 
P2 = Lfc/L^l-cxp(^-^^ (5) 

Figure n shows that the agreement with the real num- 
ber of shortcuts for dense SWGs is very good. Thus, 
in what follows, to translate the results of our model to 
dense small world graphs, it suffices to take p2 and in- 
vert Eq. Q to obtain the corresponding p. Notice that 
for values of p close to 1 there is some " overshooting" , as 
the resulting SWGs have more shortcuts than the corre- 
sponding random graphs (i. e. p2 > 7). 

III. TOPOLOGICAL PROPERTIES 

A. Average Distance 

For random graphs with a number 0{N^) of edges, it 
is known that the diameter (i. e. the largest distance 
between any two nodes) is equal to 2, for large values of 
iV. This means that from any vertex only 1 or 2 steps are 



B. Clustering coefficient 

The average distance can be thought of as a global 
property of a graph: it gives an idea of how far apart 
is any vertex from any other. The clustering coefficient 
provides a different kind of information: it is a measure 
of how locally connected is the graph. It is obtained by 
calculating, for every vertex of the graph, the number of 
links joining points of its neighborhood divided by the 
total number of possible links in the neighborhood, and 
taking its average over all vertices. This gives the prob- 
ability that two neighbors of a vertex are connected to 
each other. It can be shown that this is equivalent to 
calculating the total number of triangles in the graph, di- 
vided by the total number of paths of length 2 (hereafter 
called 2-paths). 

Instead of calculating this number for every graph and 
then average over the ensemble, we calculate an 'an- 
nealed' version of it [l^: we take the average of the num- 
ber of triangles and divide it by the average number of 
2-paths. For large values of N the agreement with the 
numerical values turns out to be very good (see Fig. [21). 

To calculate the number of 2-paths, one has to find the 
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FIG. 2: Average random distance and clustering coefficient 
for graphs Gpj(7), scaled by their maximum and minimum 
values. The symbols are averages taken on 30 matrices with 
A'^ = 6001 and 7 = 0.1. The line is the theoretical prediction 
for the average distance. Error bars are smaller than the 
symbols. 



n(f)[p?/(x) + 2pw2{\Tu\ - I{x)) 
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pUn ~2\rk\ + 1 (x))] 



(8) 



with Tl{x) = 1. The first term in the square bracket 
counts the number of paths with two short edges, the 
second counts the paths with one short and one long edge 
and the third counts those paths with two long edges. To 
calculate the number of triangles one proceeds in exactly 
the same way, but it must now be assumed that the end- 
points of the paths arc connected by an edge. Thus, the 
number of triangles is given by Eq. © using 11 (a') = l—pi 
for f e rfc(O) and U{x) = 1 -pa for f G F - rfc(O). After 
some algebra the clustering coefficient obtained is 
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with 



7dbi - P2)'((3/4)'^ - 7d) + ((pi - p2hd + P2f 



Id 



d 



(9) 



(10) 



number of 2-paths for every matrix, weight it with the 
probability of that matrix, and then sum over all possible 
graphs. An equivalent way of doing this is to sum over 
all possible paths, with a weight given by PiP2~\ where 
i is the number of long links in the path. 

In the limit of infinite N all vertices are equivalent, 
and it is enough to consider all the paths beginning from 
a fixed vertex, and then multiply by N . To count these 
paths, one needs to know the size of the intersection of 
the /c-neighborhoods of two points. If 7 < 2"^^, this in- 
tersection is a connected region and the calculation of its 
size is very simple. We analyze only this case, because 
one is usually interested in low dimensions. For bigger 
values of 7 the calculations are straightforward but much 
more cumbersome. 

Assuming the fixed vertex located at the origin, and its 
other end at a position x £ Z"^, the size of the intersection 
of the two fc-neighborhoods is 



/(f) 



UtilM-^ x€r2k{0)~Tk{0) (6) 

xer^kio) 



where TkiO) is the set of nodes that form the k- 
neighbourhood of the origin, and 



h{x) = 2ki - + 1 



(7) 



Using this, we obtain that the total number of 2-paths 
is, to first order in N, 



In the case of a graph Gp-^ (7) the clustering coefficient 
simplifies to 



C^{pi-p2f{ii/if-ld)+ld (11) 

A comparison with data obtained from graphs with 
6001 vertices can be seen in Fig 13 Theory and data only 
differ where 1 — pi ~ (jN)"^, where our approximations 
do not hold. 



IV. SPECTRUM OF THE GRAPH 

For a graph with no loops and with no multiple edges, 
as those treated in this article, the adjacency matrix A is 
a 0-1 matrix where Aij = 1 if there is an edge connecting 
vertices i and j, and otherwise. Diagonal elements are 
set to 0. The eigenvalues of A provide a lot of information 
about the structure of the graph (see For example, 

if the graph is regular and its eigenvalue distribution is 
symmetric, then the graph is bipartite. Eigenvalues also 
provide useful bounds to quantities as different as the 
diameter of a graph |l5l | and the mixing rate [T^ of ran- 
dom walks. In the following section we summarize some 
results for the spectrum of circulant graphs which will be 
useful in our determination of the spectrum of DPRGs. 



A. Circulant graphs 

The family of circulant graphs is one of the few for 
which all the eigenvalues and eigenvectors can be calcu- 
lated analytically The adjacency matrix of such a 
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graph is a circulant matrix, where row i is simply the 
first row shifted i places to the right. In our case, the 
adjacency matrices of the graphs with pi = 1 and P2 — 
are circulant. They are symmetric matrices with null di- 
agonal elements, satisfying, for d — 1 and j > i, Aij — 1 
for j — i < fc/2, and otherwise. 

The eigenvectors of this matrices are N vectors Vj = 
(l,Pj,p^2 ...,pf-^) for < J < TV - 1, with p, = 
exp (27rj/iV). Interestingly, all circulant matrices have 
the same eigenvectors. This is not the case for the eigen- 
values, which satisfy = J2f=i ^jiP)~^ ■ '^'^^ case, 
this gives, 



^ - P. 



-(fe/2+l) 



N 



,,,,^!2Ml±21f2^ti^ (12) 



N 



sm{'Kj/N) 



Except for the first one, the Perron-Frobenius eigen- 
value Xq = k, all the other eigenvalues satisfy = Xn-j 
(we consider only odd values of N). Thus, there are 
{N ~ l)/2 eigenvalues with multiphcity equal to 2. For 
j ^ N we have. 



A| ^ sin(7rj7) sin^(7rj7/2) ttj sin(7rj7) / f_ 
N TTj N 3iV2 \m 

(13) 

For j — 0{N), it can be seen from Eq. (fT^ that Xj = 
0(1). 

For d > 1 our substrates are reticles where the range 
of connection of each node depends on the spatial direc- 
tion i, through a variable ki. The adjacency matrices for 
these graphs are simply the Kronecker product (or tensor 
product) of the corresponding N^^'^ x N^^'^ matrices for 
each direction: A = A{ki) (g) A{k2) ■ • ■ A{kN). The 
resulting matrix is circulant, and has the nice property 
that its eigenvalues are 



^^n...u=I[^, for l<j.<iV (14) 
1=1 

where the A^.s are given by Eq. (|12|l . 

B. Dense partially random graphs 



by 



The moments Mj of an x matrix A are defined 



j-paths that connect vertices i and k. Thus, the trace of 
the jth power of A is the total number of closed paths, 
called cycles, of length j. 

If, for an ensemble of matrices, the calculation of the 
average number of cycles can be performed for every 
length, by using the averaged version of Eq. (|15|l it is 
possible to obtain the average distribution of eigenval- 
ues. This is the route taken by Wigner 18] in his famous 
derivation of the semicircle law for random matrices. His 
work was afterwards extended |l9j| to prove that for large 
matrices, the distribution of eigenvalues of almost ev- 
ery matrix of the ensemble tends to the semicircle. In 
Ref. [23| the same method was used to calculate the mo- 
ments of random 0-1 matrices. 

In both these works, and in many others, all the cy- 
cles are characterized as sequences of vertices. For each 
sequence of length j the statistical weight is the same, as 
all the edges in the graph have the same probability of 
being present. 

In our model this is different because short and long 
links have different probabilities. Therefore, in our case it 
is better to characterize each cycle as a succession of dis- 
tances. We begin by analizing the unidimensional case, 
and afterwards we indicate the modifications necessary 
to extend the results to higher dimensions. 

On a ring, vertices can be consecutively labeled with 
numbers from 1 to iV by going round the circle. When 
standing on a vertex, a walker can only make a step to the 
right or to the left. The distance between vertices i and 
j is defined as dij ~ min(|i — j|, TV— |i — j|). Intuitively, it 
is the shortest number of consecutive vertices, including 
the end vertex, that a walker must traverse to go from i 
to j. To define a direction we say that if a walker, using 
an existing link, goes from vertex i to vertex j such that 
(«— j) mod iV < (A^—l)/2 then it has made a step to the 
right covering a distance dij. Otherwise, we say that he 
has traveled to the left, covering a distance dij. Thus, a 
j-path can be defined by a succession of j distances. The 
total distance traveled is the sum of the signed distances. 

If the path is a closed one, i. e. a cycle, the total 
signed distance must be equal to mN, because the cy- 
cle can contain m complete rounds of the circle. In a 
cycle traversing i long edges, called a jz-cycle, the num- 
ber of complete rounds will be bounded by rumax = 
L(j - 2 + + 1)/(2^)J > where [x\ gives the largest 
integer smaller than or equal to x. Pji is defined as the 
number of ji-cycles for a fixed position in the cycle of 
the i long links. Because the probabilities for each edge 
are independent, Pji does not depend on the actual po- 
sitions of the long links. As the probability of a ji-cycle 
is p{~*p|, the average over t/p^pa (7) of the jth moment is 



N 



Mj = N-^Tr{A^) ^N-^J2 K 



(15) 



The jth power of an adjacency matrix gives informa- 
tion about all the possible paths of length j in the asso- 
ciated graph. More specifically, {A^)n: is the number of 



4=0 



(16) 



As explained in the introduction, all the quantities cal- 
culated in this article are averages. Thus, hereafter we 
drop the overlines to avoid overloading the notation. 
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Let us call di the distance traveled in the Ith step. To 
long links there correspond signed distances satisfying 
k < \di\ < N/2. For short links, the corresponding dis- 
tances satisfy 1 < \di\ < k. Thus, Pji is the number of 
solutions d € Z-* of the problem 



^di=mN for I m| < m„ 



1=1 

l<\di\<k/2 iiO<l<i 
k/2 <\di\ <{N -l)/2 \ii<l<j (17) 

Using the principle of inclusion-exclusion , Pji can 
be written as 



1=0 



(18) 



where Pji is the number of solutions to the simpler prob- 
lem 



di = mN for \m\ < 



1=1 

l<\di\<k/2 iiO<l<i 
l<\di\<{N-l)/2 ili<l<j. (19) 

Pji corresponds to the number of cycles where j — i 
fixed steps use short links, whereas the other steps can 
use either long or short links. It is not difficult to see 
that Pji is of order iV^^^ (as the path is closed, only 
j — 1 steps may be freely chosen). Notice that d; cannot 
be zero because loops are not allowed. But the presence 
of loops only adds to Pji terms of order jN^~'^ . As we are 
only interested in the dominating term, we include the 
loop terms, that allow for simpler calculations. Using 
this, and rescaling the distances, the problem of Eq. ifT^ 
can be rewritten as the number of solutions of 



where Q{x) is the Heaviside (or step) function(8(a::) = 1 
for a; > 0, and 8(x) = otherwise) and 



/Ar(p,Z,m) = J-m+p-; + 7(^-p)-^^ (22) 

As we are only interested in the dominant terms of Pij , 
we develop it in powers of TV to get 



min{l.j — i) 



m=—m,Tnax 1=0 



p=0 



p Jv-p. 



U 1) 



N- 



(^._2),(V-^)+^(^"^)j®(/-) 

(23) 



where 



foo = foo{pj,v)='--m + p-l + 7(^ - p). (24) 

For the first order term of Eq. H23|) . the sums can be 
performed (see Appendix) to obtain the simple result 



7VJ-iy-» for 1 < i < , 



(25) 



For PjQ there is not such a simple expression. On the 
other hand, Pjo = Pjo- Thus, Pjo is the number of j- 
cycles of a graph with only short links, that is, a circulant 
graph. Thus, it can be written as 



N-l 



PjO ~ Pjo - E 



(26) 



where the As are the eigenvalues of a circulant matrix, 
given by Eq. . 

Using all this, the number of ji-cycles can be calcu- 
lated, giving 



3 

Y,di= + 1) + U - i){N + l)/2 + mN 
1=1 

for I TO I < TTlmax 

l<di<k+l iiO<l<i 

l<di<N ifi<l<j. (20) 

Using again the principle of inclusion-exclusion, we get 

m^aa: j min{l,j-i) A / ■ \ 

E D-D' E cr C-> 

m=-m,„a^ 1=0 p=a \ ^ / \ -f/ 

(N fN{p,l,m) + j - l\ 

( ]Q{fN(p,l,m)), (21) 



P,, - (-l)*(P,o - N^-'j^) + N^-^j^-\l - 7)' (27) 

Using this and Eq. (|16|l . we obtain, for the first order 
of the moments of a dense partially random matrix: 

Mj = (Pjo - l'Ni-'){pi - P2y + liN^-^ for J > 3 

(28) 

where 7* = pi7 +]32(1 — 7) is the average degree of a 
vertex. The spectrum of eigenvalues that generates the 
moments corresponding to Eq. (|28|l . for all values of j 
(i. c. not only for j > 3), satisfies Xi/N 7* and 
\j/N (pi - p2)Xj/N for j > 1. Notice that this 
would imply that the whole spectrum is only a rescal- 
ing of the spectrum of the substrate (only the first eigen- 
value is scaling differently). But the problem is that. 
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FIG. 3: Average distribution of eigenvalues for Gp^ (7) with 
7 = 0.1 and pi = 0.6, for = 3001 (dashed line), N = 6001 
(dotted line) and A'^ — 9001 (full line). The averages were 
taken over 100, 30, 10 matrices respectively. The small verti- 
cal bars over the horizontal axis show the first order predic- 
tions for the discrete part of the spectrum, extended to all 
values of A. 



because of the symmetry of the matrices, the second mo- 
ment is, to first order, M2 = Nj^,, clearly different to 
what would be obtained by setting j = 2 in Eq. 
{Mj=2 = iV((7 - 7^)(pi - P2 f + 7*)). This means that 
not all the eigenvalues can tend to the values mentioned 
above. 

On the other hand, the fact that Mj = 0{N^-^) im- 
pUes (see Eq. (fT^ I that every eigenvalue that satisfies 
|A| = 0{N) must tend to (pi -_P2)A'=. And from Eq. lO 
it can be seen that the number of such eigenvalues di- 
verges with A^. But the number of eigenvalues that do 
not tend to Ac, which satisfy |A| — o(N), must also di- 
verge with A'^; otherwise, its influence would not be felt 
in M2. 

In Fig. O we show the distribution of eigenvalues for 
matrices with A^ — 6001. It can be seen that the largest 
eigenvalues are very close to the values predicted, and the 
deviations are larger for smaller values of X/N. But for 
small values of A the spectrum seems to be continuous 
and similar to a semicircle distribution. 

If we assume that the distribution is given by a discrete 
part, where to dominant order Xj ~ A^ and a continuous 
part given by a semicircle distribution (to dominant or- 
der), the second moment can be used to determine the 
width of the semicircle. It is known that the semi- 
circle distribution 

(-^Vcr2 _ ^2 foj. _ ^ < < 
(29) 
otherwise 



FIG. 4: Scaling of the average of the continuous part of the 
spectrum for Gpj^{'y) with 7 = 0.1, for different values of pi 
and N. The averages for A = 3001, A^ = 6001 and A^ = 9001 
have been taken over 100, 30 and 10 matrices, respectively. 
The full line shows the theoretical prediction (semicircle dis- 
tribution). 



generates moments 

Ml^. ^JdX X'^ps (A) = "1(1^ (^/2)'^ (30) 

(odds moment vanish because of the symmetry.) From 
Eq. (|28|) . extended to j = 2, we know that to obtain 
the correct value of the second moment for the whole 
distribution, the continuous part should satisfy = 
^(7* - 7* - (7 - 7^)(Pi - ^2)^)- For this, the width of 
the semicircle must be 



a = 2v/iV(7, - 7? - (7 - 7')bi - ^2)^) (31) 

To test the correctness of this assumption, we have 
used it to scale the average of the continuous part of 
the spectrum of Gp^ (7) for several values of A^ and pi . 
The results, displayed in Fig. ^ show that this scaling 
collapses all the curves onto the unit semicircle, as pre- 
dicted. Notice that the deviation from the semicircle is 
only noticeable for large values ofpi, because in this case 
the discrete part of the spectrum contains a significant 
part of the eigenvalues, thus partially depleting and skew- 
ing the semicircle. 

To confirm analytically the existence of the semicir- 
cle, one should develop the moments to an order large 
enough in the corrections to see the contribution of the 
semicircle. The problem is that, whereas the order of 
magnitude of AI^ is N^^^, the semicircle generates mo- 
ments Mj — 0{N^/^). Thus, for high values of j the 
contribution of the semicircle lies deeply buried under 
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the ones of the discrete part of the spectrum. Neverthe- 
less, we have calculated the moments to the second order 
(see next subsection) and confirmed that the semicircle 
provides the right value for M4. 

The existence of the semicircle imposes a limit on the 
number of eigenvalues in the discrete part: it must con- 
tain only the eigenvalues that satisfy |A| > cr. Eq. H13|) 
implies that the number of such eigenvalues is propor- 
tional to N/a = 0{Vn). 

Putting all this together the complete eigenvalue dis- 
tribution of DPRGs is, to first order in N, 



A2 



TTj I 



-a <\<a 

otherwise 
(32) 



where A^^- is the number of values of j such that 



P2)- 



3(7rj7)) 



> (7, i. e. the number of eigenvalues con- 



tained in the discrete part. It must be remarked that 
even though these eigenvalues are degenerated this de- 
generacy breaks down for finite values of N . Nevertheless 
in our simulations their separation is so small as to make 
them statistically indistinguishable inside each peak, for 
the values of TV chosen. 



1. Second order calculation 

To go beyond the dominant order in the calculation of 
the moments, two different contributions must be taken 
into account. 

The first contribution is simply a refinement of the 
calculation of Pji for i > [23, which can be split in 
two terms. One of them is simply what one gets when 
using the second order term from the development of 



i(2) 



The 



Pji (see Eq. This gives P^-' = 

other part arises when we subtract from Eq. H21|l the 
paths with loops, i. e, the solutions of Eq. (|19|l having 
di ~ for some values of I (that were introduced to 
make calculations easier). The number of such solutions 
is 0(A^^~'~^). Thus, for the second order calculation we 
need to subtract from Eq. H21|l the number of solutions 
with only one vanishing distance, which is iPj-ii-i + {j — 
i)Pj-ii. Adding these two terms, and using Eqs. H18() 
and H16() . the first contribution to the correction of the 
moments is obtained: 



am/ ^ APjo{pi-p2y -p2jMj_i 



(33) 



where APjo is the first correction to Pjq. 

The second contribution originates in the fact that the 
weight assigned to each cycle in Eq. ifTCI) . p-{~^p2, implies 
that all the edges traversed in each cycle are different. To 
correct this, some paths have to be reweighted. But the 
number of j-cycles with r edges repeated is 0{N^~^~^) 
or smaller, thus we only need to consider cycles with one 
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FIG. 5: Schematic representation of the possible cycles with 
only one repeated edge. The legends indicate the number of 
edges present in each subcycle. The edge traversed twice is 
shown as a thick line. 



edge repeated. But there are only two possible classes 
of such cycles, as shown in Fig. |5l the first class con- 
sists of two cycles sharing an edge which is traversed in 
only one direction, and the second class consists of two 
cycles joined by an edge which is traversed in both direc- 
tions. For the first class (Fig.[S^), the number of cycles 
is 0{Ni-i-'^)0{Ni-^) = 0(iVJ-3) for aU possible val- 
ues of q. Notice that in this case q must be positive for 
the edge to be traversed twice. The number of possi- 
ble j-cycles for the second class (Fig. Oa) is proportional 
to 0{N^-i-^) ■ N ■ 0{m-^) = 0{m-^) if g > 2 and 
0{N^-^)N = 0(iVJ-2) if g = 0. Therefore, the con- 
tribution of dominant order, given by those paths with 
q = (shown in fig Fig. [S]:), can be written 



AjjNMj^2 for j > 5, 



(34) 



where A7 = 7(^1 — pf) + (1 — 7)(p2 — pD- 3^ is the 
number of possibilities for the choice of the edge that 
will be repeated, the two terms in A7 correspond to the 
possibility that the repeated term is a short or a long 
link, and Mj_2 is the number of cycles of length j — 2. 

Notice that to obtain Eq. 1)34(1 we have only corrected 
the weight of the repeated edge, assuming that in the 
closed part of the paths considered (corresponding to the 
ovals in Fig. |SJ all edges are different. For j > 4 this 
is correct, because the number of such subcycles needing 
reweighting is of second order respect to the total. But 
this is not so for j = 4. In this case, the oval in Fig. 
corresponds to a single edge traversed twice. Taking this 
into account, the right contribution to the correction of 
M4 is 



AMf = 4iV^(7^-7,,), 



(35) 



where 7*, = plj + p^l — ■y). 

Adding both contributions, the first correction to Mj 
can be written as: 
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AM4 = {pi~P2f{AP4n{pi-p2)-4p2{P30-N^l^)) 

+N\2{jl-jl)-Ap2-f!) 

(36) 

Api~P2y~^jP2{Pj-io - N^-^i'-') 

+ (pi-p2)^AP,o forj>5 (37) 

Consider now the distribution of eigenvalues that has 
been proposed to generate this moments. We must in- 
troduce corrections to it, to account for the calculated 
corrections to the generated moments. For j > 4 only 
the corrections to the discrete part of the distribution 
are needed. In terms of the fist order term of the de- 
velopment of the eigenvalues, the correction to the jth 
moment is 

N 

AMj=jY,Xr^AK. (38) 

1=0 

Replacing this in the left side of Eq. l|77|l and using 
that the same holds for the corrections to the moments 
of the substrate (AP^o = AMj), a comparison of terms 
in both sides of the resulting equation gives 

AA, = AA°~p2 + AriA7, (39) 

But, by construction, these corrections to the discrete 
eigenvalues generate the right corrections only to the mo- 
ments of order larger than the fourth. For the fourth 
moment this discrepancy must be bridged by the corre- 
sponding moment generated by the continuous part of 
the spectrum. But, 

AM,=4 - AAf4 = 2(7, - 7,2)' (40) 

which is exactly the fourth moment generated by the 
semicircle distribution given in Eq H29|l . 



C. Higher dimensional substrates 

The results obtained in the preceding sections can be 
extended to DPRGs defined on higher dimensional sub- 
strates. These substrates are defined as hypercubic lat- 
tices in d dimensions where each node is connected to 
a hypercube of kik2---kd other nodes. As already men- 
tioned, we assume that the hypercube is closed, in the 
sense that nodes at the boundaries of the hypercube are 
considered nearest neighbors of the nodes at the opposite 
boundary. 



In a d dimensional DPRG a ji-cycle can be represented 
as a succession of j distances d i G Z'^. Component- 
wise, this can be regarded as the superposition of d uni- 
dimcnsional subcycles. Notice that the number of short- 
cuts used in each subcycle is smaller than or equal to i. 
To build a d-dimensional jj-cycle out of unidimensional 
paths, only the following condition must be satisfied. Let 
us call {i} the set of steps of the d-dimensional ji-cycle 
that traverse shortcuts, and {i;} its analog for the sub- 
cycle in the ^th dimension {1 < I < d). The condition is 
then that the union of all the sets {ii} be equal to {i}: 

Therefore, counting the number of possible ji-cycles 
is equivalent to counting the number of unidimensional 
subcycles satisfying this condition. This can be written 
as 

with si = I]q=i «9 and su = I]g=i ^qi- Pji is the num- 
ber of ji-paths in the Zth dimension. Notice that iiq is 
the size of the intersection of sets {ii} and {iq}- The 
evaluation of this expression is cumbersome but straight- 
forward, giving 

p,M = (-i)'(P,o(d) - iv^-^7;^) + iv^-^7r^(i - ^dY 

(42) 

where Pjo{d) = ULi ^^^^^ ^'^ = nf=i 7(0- 

Using this and Eq. Hlt)|) . we obtain for the first order 
of the moments 

Mj = {Pjo{dyi'dN'-'){Pi-P2y+{piid+P2{i-id)yN^-^ 

(43) 

for j > 3. Using the same reasoning of Section IIV Bl 
we see that the eigenvalues satisfying |A| — 0{N) must 
tend to those of the substrate, rescaled by the disorder: 
- {Pi -P2)>^'j,j,...j,, where >^,j,,„,, are given by 

Eq. IBI). 

In analogy to the one dimensional case, for the smaller 
eigenvalues we can conjecture the presence of a continu- 
ous distribution following a semicircle law, because of the 
discrepancy between the real second moment and the one 
generated by the discrete distribution. Unfortunately, 
this conjecture cannot be tested for all values of d in 
the same way used for the unidimensional case in Sec- 
tion ^^^3 The reason of this is that, as the moments 
are obtained as sums of products of d unidimensional mo- 
ments, their development involves powers of N^^'^. But 
the moments generated by the semicircle are 0{N^/'^). 
Thus, for M4 its contribution must be searched in the dth 
correction to the real fourth moment, whose evaluation, 
even though straightforward in principle, gets extremely 
cumbersome even for small values of d. 
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Fortunately, for most applications one only needs the 
largest eigenvalues in absolute value (see Sec.[V]for some 
examples), which are given by the discrete part of the 
distribution. And if a function of all the eigenvalues is 
needed, it can always be rewritten as a series involving 
the moments. 



D. Comparisons 

The limit of small 7 should give us an idea of the ap- 
proximate form of the spectrum for sparse SWGs. In this 
limit, Eqs. 10} and js} givepi ~ 1— p andp2 — 7P- From 
Eq. (|32|l . one can see that, at least for not too large values 
of j, Xj ^ piXq, for small 7. Thus, the eigenvalues accu- 
mulate at a distance of p'yN from the Frobenius-Perron 
eigenvalue Aq = 7iV, with a trail of eigenvalues reaching 
to the edges of the semicircle. In the small 7 limit, the 
width of the semicircle is cr ~ 2y/2Nj{2p — p^). 

For small values of p the eigenvalues accumulate so 
close to the Perron-Frobenius eigenvalues that the gap 
should only be visible for very large values of N . The 
continuous part of the distribution, whose width is pro- 
portional to y/p, gets very small and contains few eigen- 
values, so its shape becomes very irregular and skewed 
to the negative side (to retain the vanishing of the first 
moment). This picture is very similar to what can be 
seen in Figure 3b of Ref . 0] ■ 

If p is not small, the accumulation point is clearly sep- 
arated from the Frobenius-Perron eigenvalue. Besides, a 
can be close to \/7fy, thus including enough eigenvalues 
to take a shape close to the semicircle. This shape should 
also be skewed. This picture is very similar to what can 
be seen in Figure 3c of Ref. 3 . 

It is also similar to what was find in Ref. [^J. In 
the graphs considered in that article links are added to 
a sparse substrate (i. e. they are not rewired). The 
dense version of this corresponds to taking pi = 1 in our 
model. Even though it is the spectrum of the laplacian 
matrix that is studied, the results can be translated very 
easily to the spectrum of the adjacency matrix, because 
for large sizes the graphs can be considered regular, in 
which case the eigenvalues of both matrices can be related 
by the formula = k — X"^. It is found that two peaks 
appear. The closest to the Perron-Frobenius eigenvalue 
is separated from it by a pseudo-gap (i. e. an interval 
where there are eigenvalues, but very few of them). This 
peak is found to be "in quantitative agreement with the 
ring spectrum" , and can be related to the accumulation 
point mentioned above. The other peak, which is very 
irregular for small number of shortcuts, can be related to 
the continuous part of the spectrum found in DPRGs. 

V. SOME APPLICATIONS 

It is interesting to notice that the average distance and 
the clustering coefficient present qualitatively the same 



behavior as that seen in sparse SWGs. We can see that 
for small values of pi the graphs obtained have aver- 
age distances which are close to those in random graphs, 
while retaining a clustering coefficient close to the val- 
ues present in circulant graphs. Naturally the range of 
values spanned by (the logarithm of) both quantities is 
much larger in SWGs. 

Having the distribution of eigenvalues, or, equivalently, 
the expression for all the moments, allows one to calcu- 
late, or at least to bound many processes that can take 
place in DPRGs. 

For regular graphs (i. e. graphs where all vertices have 
the same degree), the spectrum of the adjacency matrix 
can be very simply related to the spectrum of the lapla- 
cian matrix, defined hy L = D — A where A is the adja- 
cency matrix and D is the degree matrix (a diagonal ma- 
trix such that da is the degree of vertex i) and the normal 
matrix N = D~^A. The laplacian has very interesting 
properties and many applications in physics, specially 
because it arises in the discretization of the Laplacian 
operator As SWGs are regular, their laplacian and 
adjacency matrix eigenvalues are related by = k — X^. 
If the eigenvalues are ordered from small to large, the first 
eigenvalue is Aq = 0. The second eigenvalue is probably 
the most important as it can be related to a number of 
properties of processes taking place in such graphs. 

In the following we show a couple of examples where 
we apply the results obtained in the preceding sections. 

A. Mixing rate 

A random walk on a graph is defined as a Markov chain 
where the probability of jumping from vertex i to vertex 
j is 1/dj if they are connected, and otherwise 
eral properties of random walks can be related to the 
spectrum of a graph. For every time t there will be a dif- 
ferent probability Pt{j) of finding the walker on a site j. 
A stationary distribution is defined by the requirement 
that Pt = Pt+i for t > T, for some T. 7r(j) = dj/2M 
is a stationary distribution, and for regular graphs it is 
unique. It can be shown that, regardless of their initial 
state, all the walks tend to this distribution, provided the 
graph is connected and not bipartite. 

When the walk has reached the stationary distribution 
it has essentially lost all memory of its initial state (or 
distribution) and all the vertices are sampled with prob- 
ability proportional to their connectivity, which is useful 
for several algorithms. But how fast is the convergence 
to the stationary distribution? One of the possible mea- 
sures of this is the mixing rate, defined as 

fi = limsupmax \p,j - dj/M\^/* (44) 

It can be shown that /i — A^, the largest nontrivial 
eigenvalue of the normal matrix. Using Eqs. (|32|l and Q 
we obtain. 
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(pl - P2) 



sin(7r7) 



7r7 

7 — 1 + exp 



7P \ \ sin (777) 



1-7 



7r7^ 



(45) 



where the last equahty is vahd for a dense SWG. This 
shows that for fixed and small 7 the mixing rate decreases 
almost linearly with the disorder. Notice that the fact 
that the average distance jumps to its minimal value at 
p = 0+ docs not influence the mixing rate. 



B. Synchronization of coupled oscillators 

One of the most interesting processes that can take 
place on a network is the collective dynamics of an array 
of coupled oscillators. And perhaps the most striking col- 
lective state is that where all the identical oscillators get 
synchronized. Naturally, synchronization is not always 
possible, it depends on the specific properties of the os- 
cillators as well as on the topology of the network. In 
Ref. |23j| a very useful formalism was introduced to study 
the conditions for the existence of a stable synchronized 
phase, for a wide class of oscillators and couplings. The 
equations of motion for the ith oscillator in the network 
are 



F(x, 



N 



H(x,) 



(46) 



where F governs the dynamics of each individual oscil- 
lator, H is an arbitrary output function, a gives the 
strength of the coupling, and L is the laplacian matrix 
of the network. It can be shown that for a system 
of this form, the condition for the existence of a stable 
synchronous state reduces to 
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FIG. 6: First ten eigenvalues of the discrete spectrum of 
the adjacency matrix of a DPRG. The dashed line (A = 
~2'yN{pi — p2)/37r), that joins the first minima of all the 
eigenvalues, is a reasonably accurate estimation of the small- 
est eigenvalue for every 7. 




0.20 



FIG. 7: Synchronization threshold for graphs Gpi(7). The 
inset shows the critical fraction of links that must be short- 
cuts. 



Ai/Ai < a (47) 

where A£ and Ag are, respectively, the largest and the 
smallest nontrivial eigenvalues of the laplacian. /? is a 
parameter that depends only on the oscillators and its 
coupling, and not on the topology. (5 G [5, 100] for several 
chaotic oscillators E^. 

Using Eq. (|47|l we can calculate the synchronizabil- 
ity threshold for dense SWGs. For DPRGs and suffi- 
ciently large values of N the smallest and largest (ex- 
cluding the Frobenius-Perron) eigenvalues of the adja- 
cency matrix are located in the discrete part of the 
spectrum. The largest nontrivial eigenvalue is always 
\^ = N{{pi —P2) sin(7r7)/7r, but the index of the smallest 
eigenvalue depends on 7 (see Fig.|S|) An approximation 
to it is given by Xg = —N^{pi — p2)2/37r, which is accu- 
rate enough for our illustrative purposes. 



The synchronizability threshold is defined as the small- 
est value of p for which the system becomes synchro- 
nizable. Using the already mentioned relation between 
the adjacency and laplacian matrices, the synchroniza- 
tion condition for DPRGs becomes 

/^-l , 2 ^ /?sin(77r) ^^g^ 
Pl — P2 Svr 77r 

Notice that dense random graphs are always synchro- 
nizable, as Xl/Xs — 1. As was done in Ref. 
we have calculated the synchronization threshold for a 
dense SWG for the case of a system of oscillators with 
/3 = 37.85. The results are displayed in Fig. [71 

The region below the curves gives the set of param- 
eters for which the system is not synchronizable. It is 
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interesting to notice that this impHes that to obtain a 
synchronizable system it is not enough to have a macro- 
scopic number of shortcuts. Also, the fact that in this 
region p is positive shows that the onset of synchroniza- 
tion does not depend on the average distance because, 
as we have seen in Section Fill Al for p > the average 
distance is the smahest possible. This supports the idea 
that average distance alone is not a relevant factor for 
the onset of synchronization. But in Ref. |25[ it has been 
argued that synchronization seems to depend on the com- 
bined effects of small average distance and uniformity of 
the connectivity distribution. Yet, our results show that 
it must at least depend on some additional factors. For 
fixed values of 7 we have seen that the synchronization 
threshold is > 0, even though all graphs Gp{j) have 
the same connectivity {Nj connections per node) and 
the same average distance for p > 0. 



of the peaks can in principle be associated to large eigen- 
values of the substrate, all rescaled by the same value 
(dependent on the disorder), and the other to a continu- 
ous distribution of small width that is usually a signature 
of disorder. 

We have applied our results to the calculation of the 
mixing rate of a random walk on the graph, and to the 
calculation of the synchronization threshold of a system 
of coupled oscillators placed on the nodes of the graph. 
We have shown that below and above the threshold there 
exist graphs with the same average distance and the same 
connectivity (i. e. the same number of connections per 
vertex). Previously it has been argued that a combina- 
tion of these two factors was responsible for the onset of 
synchronizability. Our results show that there must at 
least exist more quantities involved. 



VI. CONCLUSIONS 

In this article we have studied some properties of dense 
partially random graphs, that can be obtained by adding 
edges with probability p2 to a dense ordered substrate, 
and then deleting the original edges of the substate with 
probability l—pi- We have found that they show qualita- 
tively the same behaviour as the corresponding properties 
of sparse graphs, in particular small world graphs. For 
the average distance we find that, for any macroscopic 
number of shortcuts, it falls to its minimum possible 
value (in the infinite size limit.) We show that the clus- 
tering coefficient decays slowly, thus allowing for a range 
of parameters where the graphs have relatively large clus- 
tering and minimal average distance. 

By counting cycles on the graph, we have obtained the 
distribution of eigenvalues of the adjacency matrix. We 
found that it consists of two parts: a discrete one where 
the eigenvalues, of order N , are simply rescalings of the 
corresponding eigenvalues of the substrate, and a contin- 
uous part given by a semicircle distribution whose width 
is of order \/N, and is proportional to the disorder. It is 
interesting to notice that a similar form of the spectrum 
has been obtained |2fl] for matrices that are the sum of 
a stochastic matrix and a block matrix. The block ma- 
trix is composed of k blocks of size proportional to N 
where all the off diagonal components are equal. The 
spectrum obtained consists of a semicircle distribution of 
width proportional to ^/N and a discrete part containing 
the k largest eigenvalues (of order N). Notice that in 
our case the discrete part contains a diverging number of 
eigenvalues. 

We have shown how the distribution found can be 
useful to understand the distributions that arise in the 
numerical study of the spectrum of sparse small world 
graphs. In the studies published so far two peaks ap- 
pear, and a pseudo-gap separates the bulk of the spec- 
trum from the Frobenius Perron eigenvalue. Comparing 
with our results for small values of the connectivity, one 
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Appendix 

Here we sketch the calculation that leads to Eq. H25|l . 
Although we suspect that there must be a shorter path 
to Eq. (|25() . we have not been able to find it. 

We begin by splitting the dominant term in Eq. 12,'-{|l . 

pjP , in two parts: 



J min(/ j'— i) , 



7n— — mr,iax l—O 



p=0 



p J \l-p^ 



r(A + B) 



(49) 



where g{x) = x^^^O{x) and 



•ff(z/2-m + p-; + 7(^-p)) (50) 

m = -mrrtax 1=0 P=0 ^ ^ ^ ^ ^' 

.g{t/2-m + p-l + j{^-j^-p)) (51) 

Notice that, to avoid overloading the notation, we 
extend the definition of the combinatorial numbers, 
C) = Iz4tw b<a,to il) = for b > a. By 
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making the replacements p j — i — p and I j ~ i — I, 
and rearranging the sums, we obtain: 



m = — nimaa; P 1 = 

■g{i ~ V — p — I + 7p) = A' - B 



(52) 



where / = X!^^o(~^)*(i) operator and 



j—i i 2m„ 



A' = (-l^E'E' E 5(^^ + ? + 7P-a) (53) 

p I m— 



B' = (-IF"-' E E(-i) 

m— p— 1 

•5(1; + ^ + 7P — a) 



J - « 



^ ' 1=0 



(54) 



By manipulating the indices and reordering the sums, 
it is not difficult to show that B = B'. Using the def- 
inition of a, we can show that [I + jp ^ a\ < 2mmax- 
Thus 



A' = (-iF'E'E' E i-m+ip-ay-^ 
= E'E(-i^'O) E(-™+7P-«r' + 

p 1=1 ^ ^ m=-l 

+ E (-^r-i'-')±' 

hp-a] 

■ {—m + jp—ay'^ 

m=0 



(55) 



It is known QJ^ that 



(-l)''a!S'(n,a) for < a < n 
^/fc".= <' (56) 
k I for a > n 



where S{n, a) are the Stirling numbers of the second 
kind j23|. Using Eq. 1)56(1 one sees that the term in the 
last line of Eq. ((55|l vanishes (but notice that this can 
only happen for i > 0). Thus, expanding in powers of m, 
we get 



1=1 



A' = E'E(-i)"n E(7P-"r^-'=E 



A:=0 



m— 1 



(57) 



Bernoulli's expression for a sum of pow- 
ers 28] ^is^^^jI]L=i = Y.ttlhrnl'" where 
(^t^)Bk-m+i and Bi are the Bernoulli 



bk 



(-1) 



fc+l V k 

numbers. Using this, and summing over I, we get 



A' = E P , ^)E^"^(-1)""-''*!^(«'*) 



fc=i-l 



■ E '(^^^ " 



(58) 



But Eq. 156|) shows that in Eq. (|58l) the sum over p 
vanishes for k > i — 1. Thus, only one term survives, the 
one with fc = i — 1. Using that Bq — 1 and S{i,i) — i\ 
we finally obtain 



(j-1) 



■ A' = N^-^^^-' 



(59) 
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